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The masses of the S'-wave mesons consisting of the light {u, d, s) quarks are calcu¬ 
lated within the constituent quark model. The relativistic Schrodinger-like equation 
with a confining potential is numerically solved for the complete relativistic qq po¬ 
tential including both spin-independent and spin-dependent terms. The obtained 
masses of the ground-state vr, p, K, K* and (j) mesons and their first radial excita¬ 
tions are in reasonably good overall agreement with experimental data. 
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The theoretical description of the properties of light mesons such as vr, p, K, K* and 0 
within the constituent quark model presents additional difficulties compared to heavy-light 
mesons and heavy quarkonia. In fact, due to the highly relativistic dynamics of light quarks, 
the v/c and l/rrig expansions are completely inapplicable in the case of light mesons and the 
QCD coupling constant a* at the related scale p is rather large. Moreover, the behaviour 
of as{p^) is unknown in the infrared region and is thus model dependent (exhibiting, e.g., 
freezing behaviour, etc.). The pseudoscalar mesons vr and K produce a special problem, 
since their small masses originate from their Goldstone nature. 

Different attempts to study light mesons on the basis of the relativized quark model 
the Dyson-Schwinger and Bethe-Salpeter equations gi, chiral quark models with spon¬ 
taneous symmetry breaking (e.g. the Nambu-Jona-Lasinio model) g and lattice QGD g 
were undertaken. Here we consider the possibility of describing light mesons on the basis 
of the three-dimensional relativistic wave equation with a confining potential. While de¬ 
scribing the properties of heavy-light mesons g, we treated the light quarks in a completely 
relativistic way. Now we apply this approach to calculating the masses of light mesons in 
the framework of the previously developed relativistic quark model. 

In this approach a meson is described by the wave function of the bound quark-antiquark 
state, which satisfies the quasipotential equation of the Schrodinger type g 
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where the relativistic reduced mass is 
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Here M = £'1 + £'2 is the meson mass, mi ^2 are the quark masses, and p is their relative 
momentum. In the center-of-mass system the relative momentum squared on mass shell 
reads 
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The kernel V (p, q; M) in Eq. is the quasipotential operator of the quark-antiquark in¬ 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in¬ 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear conhning potentials, 
where the vector conhning potential contains the Pauli interaction. The quasipotential is 
then dehned by 

I/(p,q;M) = Ui{p)u 2 {-p)V{p, q; M)ui{q)u 2 {-q), (5) 


with 

v(p. q; M) = jO.B,„(k)7!‘72“ + V'l,(k)r5T2,,. + ri„,(k). 
where as is the QCD coupling constant, is the gluon propagator in the Coulomb gauge 
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and k = p — q; 7 ^ and u{p) are the Dirac matrices and spinors 
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with e{p) = \/p^ + m‘^. The effective long-range vector vertex is given by 


r^(k) = 7 ^ +( 8 ) 

where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar conhning potentials in the nonrelativistic limit reduce 
to 
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reproducing 
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where e is the mixing coefficient. 

The light constituent quark masses = 0.33 GeV, = 0.5 GeV and the 

parameters of the linear potential A = 0.18 GeV^ and B = —0.3 GeV have the usual values 
of quark models. The value of the mixing coefficient of vector and scalar conhning potentials 
e = —1 has been determined from the consideration of charmonium radiative decays [^. 
Finally, the universal Pauli interaction constant k = —1 has been hxed from the analysis of 
the hne splitting of heavy quarkonia ^Pj- states [^. In the literature it is widely discussed 
the’t Hooft-like interaction between quarks induced by instantons [^. This interaction can 
be partly described by introducing the quark anomalous chromomagnetic moment having 
an approximate value k = —0.744 (Diakonov j^). This value is of the same sign and order 
of magnitude as the Pauli constant k = — 1 in our model. Thus the Pauli term incorporates 
at least part of the instanton contribution to the qq interaction.^ 

The quasipotential © can be used for arbitrary quark masses. The substitution of the 
Dirac spinors o into © results in an extremely nonlocal potential in the conhguration 
space. Glearly, it is very hard to deal with such potentials without any additional transfor¬ 
mations. In oder to simplify the relativistic qq potential, we make the following replacement 
in the Dirac spinors: 

G,2(p) = \/"^l,2 + P^ ^^1,2 (11) 

(see the discussion of this point in j^). This substitution makes the Fourier transformation 

of the potential © local. We also limit our consideration only to the F-wave states, which 

further simplihes our analysis, since all terms proportional to vanish as well as the spin- 

orbit ones. Thus we neglect the mixing of states with different values of L. Galculating 

V s 

the potential, we keep only operators quadratic in the momentum acting on Vcoui, Konf 
replace —>■ 2 ~ 2 higher order operators in accord with Eq. m preserving the 

symmetry under the (1 <->■ 2) exchange. 

The substitution m works well for the conhning part of the potential. However, it 
leads to a hctitious singularity (5^(r) at the origin arising from the one-gluon exchange 
part (AVcouK^")), which is absent in the initial potential. Note that this singularity is not 
important if it is treated perturbatively. Since we are not using the expansion in u/c and are 
solving the quasipotential equation with the complete relativistic potential, an additional 
analysis is required. Such singular contributions emerge from the following terms 


1^2 ^ 

if we simply replace 61^2 —^ Ei^ 2 - However, the Fourier transforms of expressions m are 
less singular at r —> 0. To avoid such hctitious singularities we note that if the binding 
ehects are taken into account, it is necessary to replace 61,2 —>■ Ei ^2 — hi, 2 V, where V is 
the quark interaction potential and r]i ^2 = '^ 2 , 1 /+ ’^ 2 )- At small distances r —>■ 0 , the 
Goulomb singularity in V dominates and gives the correct asymptotic behaviour. Therefore, 
we replace ei ,2 —>■ Ei ^2 — hi, 2 hcoui in the Fourier transforms of terms (fT^ (cf. 0 ). We used 
the similar regularization of singularities in the analysis of heavy-light meson spectra j^. 


^ As is well-known, the instanton-induced’t Hooft interaction term breaks the axial Lhi(l)-symmetry, the 
violation of which is needed for describing the r] — q' mass splitting. We do not consider this issue here. 
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Finally, we ignore the annihilation terms in the quark potential since they contribute only 
in the isoscalar channels and are suppressed in the ss vector channel. 

The resulting qq potential then reads 


V{r) = Vsi{r) + VsB{r), 

where the spin-independent potential for S'-states (L^ = 0) has the form 
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and the spin-dependent potential is given by 
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Here we put Ug = aJulo) with /*i 2 = 2mim2/{mi + m 2 ). We adopt for as{fi^) the simplest 
model with freezing |l0|, namely 
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where the background mass is Mb = 2.24:i/a = 0.95 GeV Q , and A = 413 MeV was hxed 
from htting the p mass. We put the number of flavours rif = 2 for vr, p, K, K* and rif = 3 
for (f). As a result we obtain as{pl^) = 0.730, as(/*^s) = 0.711 and Q!s(/iL) = 0.731. 
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TABLE I: Masses of light S'-wave mesons (in MeV) 


Meson 

State 

n25+iLj 


Theory 



Experiment 
PDG [nj 

this work 

Ref.[lJ 

Ref. [2] 

Ref. [3] 

vr 


154 

150 

138 

140 

139.57 

P 

l^^i 

776* 

770 

742 

785 

775.8(5) 

it ' 

2^51 

1292 

1300 


1331 

1300(100) 

p' 

2^51 

1486 

1450 


1420 

1465(25) 

K 


482 

470 

497 

506 

493.677(16) 

K* 


897 

900 

936 

890 

891.66(26) 

K' 

2^51 

1538 

1450 


1470 

1460? 

K*' 

23^1 

1675 

1580 


1550 

1717(27) 

0 


1038 

1020 

1072 

990 

1019.46(2) 

0' 

2^51 

1698 

1690 


1472 

1680(20) 


* fitted value 



FIG. 1: The reduced radial wave functions R{r): (a) vr (solid line) and p (dashed line) mesons; 
(b) K (solid line) and K* (dashed line) mesons. 


The quasipotential equation m is solved numerically for the complete relativistic poten¬ 
tial m which depends on the meson mass in a complicated highly nonlinear way. The 
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obtained meson masses are presented in Tabled] in comparison with experimental data 
and other theoretical results QQi- This comparison exhibits a reasonably good overall 
agreement of our predictions with experimental mass values. We consider such agreement 
to be quite successful, since in evaluating the meson masses we had at our disposal only one 
adjustable parameter A, which was hxed from £t_B^ the p meson mass. All other parame¬ 


ters are kept the same as in our previous papers p, l7|. The obtained wave functions of light 
mesons can be used further for describing their properties (charge radii, decay constants 
etc.) and various hadronic matrix elements. The wave functions of vr, p and K, K* mesons 
are shown in Figs, ^a) and m, respectively. The 0 meson wave function looks similar. 
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